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Abstract 



The physical and bio-mechanical principles that govern undulatory move- 
ment on wet surfaces have important applications in physiology, physics, 
and engineering. The nematode C. elegans, with its highly stereotypical 
and functionally distinct sinusoidal locomotory gaits, is an excellent sys- 
tem in which to dissect these properties. Measurements of the main forces 
governing the C. elegans crawling gait on lubricated surfaces have been 
scarce, primarily due to difficulties in estimating the physical features at the 
nematode-gel interface. Using kinematic data and a hydrodynamic model 
based on lubrication theory, we calculate both the surface drag forces and the 
nematode's bending force while crawling on the surface of agar gels. We find 
that the normal and tangential surface drag force coefficients during crawl- 
ing are approximately 220 and 22, respectively, and the drag coefficient ratio 
is approximately 10. During crawling, the calculated internal bending force 
is time-periodic and spatially complex, exhibiting a phase lag with respect 
to the nematode's body bending curvature. This phase lag is largely due to 
viscous drag forces, which are higher during crawling as compared to swim- 
ming in an aqueous buffer solution. The spatial patterns of bending force 
generated during either swimming or crawling correlate well with previously 
described gait-specific features of calcium signals in muscle. Further, our 
analysis indicates that changes in the motility gait of C. elegans is most 
likely due to the nematode's adaptive response to environments character- 
ized by different drag coefficient ratios. Key words: undulatory locomotion; 
lubrication; soft matter 



Introduction 



Many undulating organisms including earthworms and nematodes move in 
environments where surfaces are lubricated by a thin liquid film. Exam ple s 
include soil dwelling nematodes that are responsible for soil aeration (|l2l ) 
and parasitic nematodes that move inside the small intestine (l25l ). The ne- 
matode and model organism Caenorhabditis elegans is of particular interest 
because it is widely used in biomedical research and usually investigated 
while crawling on wet agar plates (0, S, [20I ) . Undulatory locomotion on wet 
surfaces bears important differences with motion on dry surfaces where the 
texture, structure, and patterns of the organism's outer bod y a nd surface 
play an important role in overcoming dry friction 

(i, S HQ. On wet 

surfaces, the thin liquid film formed between the organism's outer body and 
the surface is known to reduce (dry) friction, and hence facilitates propulsion 
by providing thrust. However, many questions remain regarding the forces 
nematodes need to generate in order to overcome surface drag. Knowledge 
of such forces is critical for developing accurate dynamic models 13) , 
and understanding the nematode's internal muscle activities ([2^,0? ). It 
can also provide insight into the different kinematics and motility gaits (e.g. 
swimming and crawling) observed in C. elegans studies (0, @, 24 , IWt ) . 

Much effort has been devoted to understand the physical mechanisms 
by which organisms such as the nematode C. elegans generate thrust and 
propulsion on lubricated surfaces duj, |20|, l2J, 1271 ) . Of particular interest is 
the drag force, and its components, acting on the surface of the nematode. 
Values of surface drag forces are often obtained using resistive force theory 
(RFT) flHB, born out of studies on undulatory swimming of organisms 
immersed in a liquid (15, 1^, 2^). In RFT, the components of the drag force 
normal (-F„) and tangential {Ft) to the nematode's body are assumed to be 
linearly proportional to the local body speed along normal and tangential 
directions, respectively. The corresponding drag force coefficients are C„ 
and Ct- Under this scenario, the anisotropy between the normal and tan- 
gential drag coefficients (i.e. C„ / Ct) leads to net motion (22). For C. 
elegans swimming in fluids, the drag coefficient ratio Cn/Ct is typically in 
the range of 1.4 to 2 aaa. m contrast, estimated values of for 
nematodes crawling on lubricated surfaces (agar gels) vary by as much as 



an order of magnitude l20|, \24, l27l ). and can range from 1.5 to 35. As 



a consequence, it is difficult to accurately estimate the forces and energy 
necessary for nematodes to move on wet surfaces. 

In this study, we investigate the mechanics of undulatory locomotion on 
thin liquid films in experiments and using a simple model. By incorporating 
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kinematic data into a lubrication model, we are able to estimate values of 
the surface drag forces and the corresponding drag coefficients of C. elegans 
crawling on a wet agar surface. In addition, the nematode's internal bending 
forces during the crawling gait are estimated by considering the contribu- 
tions of internal elastic and external viscous (surface) drag forces. Finally, 
spatial-temporal signals of both kinematics and internal bending forces for 
crawling nematodes are compared to those for nematodes immersed in a 
fluid in order to gain further insight into the mechanisms governing these 
two motility gaits. The analysis provided here indicates that the adoption of 
a crawling gait appears to be a result of the specific anisotropy of the drag 
forces characterizing the wet surface rather than the high external viscous 
load. Results show that the spatial and temporal patterns of the internal 
bending force correlate well with previously described gait-specific features 
of calcium signals in muscle (0). 



Results 

Experimental observations, film lubrication and anisotropic 
drag forces 

C. elegans are transferred and allowed to move freely on the surface of Nema- 
tode Growth Medium (NGM) agar plates, which are cultured and prepared 
using standard methods (see Appendix) . After solidifying, a thin liquid film 
typically forms on top of the agar gel since the NGM water concentration 
is relatively high (98% by weight). The crawling motion of the nematodes 
is captured using video microscopy at 15 frames per second (see SVl). A 
snapshot of a single nematode crawling on an agar plate over three bending 
(undulatory) cycles is shown in Fig. 1(a), in which the red line corresponds 
to the nematode's body centerline ("skeleton"). Here, the nematode's lo- 
comotion is decomposed into normal and tangential directions; the normal 
direction is defined to be perpendicular to the nematode body length while 
the tangential direction is defined along the nematode body length (Fig. 
la). The propulsive and drag forces associated with nematode crawling on 
thin liquid films are analyzed along these two directions. 

Since the nematode's body length (L w 1mm) is approximately 25 times 
its body radius {Rw ~ 40//m), we consider C. elegans to be a long, slender 
body. Hence, the analysis (to be discussed below) is restricted to the ne- 
matode's two-dimensional cross-section (Fig. [TJj). As shown in Fig. 1(a), a 



groove forms as the nematode moves on the agar surface (Illl.l34l). In the nor 



mal direction, the nematode slides on the surface against the groove, which 
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is lubricated by a thin liquid film. The shape of this groove is described by 
a symmetric circular arc curve of radius Ra with an angular span of 29c, as 
shown schematically in Fig. [D^c). A typical groove is 1.0 fim deep and 28 
/im wide, as measured by optical interferometer (jH). The groove effective 
radius and angular span are calculated to be Ra = 100/im and 9c = 20°. The 
height of the liquid film h is defined as the distance between the nematode 
body and agar surface (Fig. [lb). For nematodes crawling on wet agar gels, 
this liquid film has been shown to be very thin such that h <C Rw (34). 



Analytical expressions for drag forces during crawling 

In the normal direction, the nematode's body and agar surface (i.e. groove) 
are approximated by two parabolic curves ha = /2Ra and = /iq + 
/2Rw, respectively (Fig. Wp)- The liquid film thickness is h = — ha = 
ho + (l/2)(i?~^ — Ra^)x'^, where /iq is the minimum film thickness. We 
define a groove shape parameter and a normalized film thickness as C = 
{1/2){R~^ — Ra^)R^ and e = h^/C, respectively. The liquid film thickness 
can then be expressed as /i = C(e + X^), where X = x/Ryj. Here, the 
groove spans from the left edge x = —Xm to the right edge x = Xm, with 
a width of 2xm = 28 fim (Fig. Ic). We note that the local curvature of 
the groove does not play a significant role on locomotion dynamics (fH. 33). 



Rather, the film thickness h determines the hydrodynamics of lubrication. 
In other words, within the lubrication region |x| ~ \/2i?^/i, the fiow can be 
effectively described by the nematode body y = h = {h^ — ha) sliding on a 
flat surface (y = 0). 

In the experiments, the nematode's typical normal sliding speed Vn is 
1 mm/s, and the viscosity and density of the liquid (water) are ^ = 1 
mPa-s and p = 10^ kg/m^, respectively. The Reynolds number, defined as 
Re = 2pVn\/2Rahl [I is less than 10~^, and the dimensionless gapjwidth 5 = 
hjRa ~ 0.01 for a typical lubrication film thickness of /i ~ 1/um (]14)- Hence, 
we may omit terms in 0{S) and 0{bRe) in the dimensionless momentum 
equations (see detailed derivation in (jsil )). If we further assume that liquid 
flow lubricating the nematode's body and the gel are uniform and steady, 
the governing equations reduce to (331 ): 

I ^ »■ 

dx ^ dy"^ ' 
duy ^ dux Q 
dy dx 
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The no-slip boundary conditions on the nematode's body and the agar sur- 
face are respectively (Fig. [TJj), 



Ux{y = 0) = 0, Ux{y = h) = -Vn cos 6 ^ -Vn, (4) 
Uy{y = 0) = 0, Uy{y = h) = -Vnsm9 = -Vn^- (5) 

By solving Eqs. 1-3 along with the boundary conditions, we obtain the 
the following velocity profiles for and Uy, 

U = —^y^ - I (v +—^\y (6) 



1 d p 3 1 d / h dp\ 2 Vn dh 



6fi dx'^'^ ^ 2 dx \2fx dx ) ^ 2h'^dx^' 

Next, we introduce the above equations for the velocity components Ux and 
Uy into Eq. 3, and obtain: 

d f h? dp\ Vn dh 



dx \ 12/i dx J 2 dx 



Equation [8] is the well known Reynolds equation ([ij, l33l ) . Note that the 
presence of a liquid meniscus has been observed on both sides of the ne- 
matode body as it crawls on agar gels. This liquid meniscus imparts an 
additional force on the nematode, namely capillary force (surface tension), 
in the direction of gravity. This capillary force enter in the analysis as a 
boundary condition (see below). 

The pressure p and its gradient dp/dx can now be computed by rearrang- 
ing and integrating Eq. El The boundary conditions are (i) dp/dx = at 
6 = 9c{x = Xm, X = Xm) and (ii) p = po+Pm at 9 = 9c {x = Xm, X = Xm), 
where p = po = is the reference pressure, and pm is the capillary pressure 
resulting from the surface tension of liquid meniscus. The capillary pres- 
sure Pm{^ —125 Pa) is estimated using the Laplace- Young equation and 
geometrical considerations (|33l). 

This results in 



dp ^ GfiVn X^ - X^ 
dx C72 (e + X2)3' 

P = -^^{Me,X)-XlMe,X)) + D, (10) 



5 



where, 

Ji{e,X)= t^/{e + t'^fdt,J2{e,X)= l/{e + t'^fdt, 
Jo Jo 

D = {6^iVnR^)/C\Jl{e,X^) - XiJ2{e,X^)) 

Given the pressure gradient dp/dx, we compute the velocity profiles Ux 
and Uy in Eqs. E] and [71 

The drag forces acting on the nematode are obtained by integrating the 
fluid (water) viscous stress (r) over the nematode's body surface, given by 
y = h = [hyj — ha). Note that r = where D is the rate of deformation 

tensor and ^ is the fluid visosity. Then, the drag forces in the normal and 
tangential directions are given by (see derivation in (j33l )) 



Fx{y = h) = fiv„ 



:Xr, 



~c 

-1 



arctan 



■ arctan 



X„ 



Xra 

7^ 



-2 + -Xle 2 



and 



Fz{y = h) = fivt 



2 arctan 



sin 9 c 



=-1/2 



(11) 



(12) 



where vt corresponds to the nematode's tangential (crawling) speed, and 
and Fz are forces per unit length. 

We note that the drag force relations obtained above (Eqs. [TT] and 
T2|) resemble the form of expressions obtained using resistive-force theory 
(RFT) (figl ). Namely, the components of the drag force are linearly propor- 
tional to the fluid viscosity (^) and the nematode's crawling speed (vn and 
Vt), such that Fx^z = tJ'Vn,tCn,t, where Cn,t is the drag coefficient and n and 
t correspond to the normal and tangential directions, respectively (see Fig. 
la). In Eqs. [TT] and [T2t the drag coefficients are defined by the expressions 
in the curly brackets. 

The force relations derived above (Eqs.[TT]and ll2p can provide some use- 
ful insight. We note that as the groove surface becomes fiatter, i.e. i?a — )■ oo, 
the shape parameter C = (l/2)(i?~^ — R~^)R^ increases, and the effective 
values of C„ and Ct become smaller. That is, the mere presence of a groove 
increases the drag forces in both normal (x) and tangential [z) directions. 
Also, the drag coefficient ratio Cn/Ct remains independent of the groove 
shape parameter C. Therefore, according to this analysis, the presence of 
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the groove is not a necessary condition for locomotion on agar plates, as 
recently suggested Moreover, any change in agar stiffness affects the 
parameter C, yet it does not alter the drag coefficient ratio, in agreement 
with previous observations (|2ol ). 



Estimating the drag coefficients 

The next step in the analysis is to calculate the drag coefficient ratio C„/Ct 
and its individual components for crawling motion. We begin by noting 
that, in resistive force theory, the ratio C„/Cj is related to the nematode's 
motility kinematics such that its normalized forward speed U/c is given 



by 4101, 120, l3i 



_ {2^^A^/\^){CJCt-l) 
""/^ {2^^Ay\'^){Cn/Ct) + l ^ ' 

The quantities c, A, and A correspond to the nematode's wave speed, bend- 
ing amplitude and wavelength, respectively. These kinematic quantities are 
measured in experiments by tracking the nematode's motion on agar plates 
as shown in Fig. 2. The nematode's motion and its kinematics are in- 



vestigated using an in- house code (j32l ). The body centroid and the body 



centerline ("skeleton") are tracked and presented in Fig. E^a) and Fig. ^h) 
in the laboratory and nematode's frames, respectively. One can extract rel- 
evant kinematic metrics from the body centroids and centerlines such as 
forward locomotion speed U and body curvature k (Fig. [2^, c). 

The curvature k is defined as k{s, t) = d(f)/ds, where (j) is the angle made 
by the tangent to the x-axis in the laboratory frame at each point along 
the body centerline, and s is the arc length coordinate spanning the head 
of the nematode (s = 0) to its tail {s = L). Body curvature is computed 
from the measured nematode's body centerlines y{s,t), where the curvature 
K{s,t) = d'^y/dx^ and s ~ 2; for small amplitudes. Fig. Wis) shows the 
experimentally measured spatio-temporal evolution of the nematode's body 
curvature K{s,t) for three beating cycles. The periodic diagonal patterns 
seen in Fig. EJ^b) show undulatory bending waves propagating from head 
to tail. The undulating waveforms are characterized by the interval of the 
periodic patterns and the slope of the diagonal lines. The body bending 
frequency (f) is obtained from the one-dimensional fast-Fourier transform 
of the curvature field k at multiple body positions s/L (Fig. 2d). A single 
frequency peak / ~ 0.5 Hz is found in the Fourier spectrum. This single 
peak is irrespective of body position and corresponds to a wave speed c = 
0.39 mm/s. The wavelength of the undulatory wave is defined as A = c// = 
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0.79 mm. A summary of the experimental results is presented in Table 1 
along with a comparison to data for nematodes swimming in buffer solution. 

By inserting the experimentally measured values of U , A, A and c into 
Eq. [T3l we estimate the value of the drag coefficient ratio for crawling 
nematodes on wet agar surfaces to be approximately C„/Ct = 9.4 it 0.6. 
This value is much larger than for nematodes immersed in Newtonian fluids 



{Cn/Ct ~ 2) (j6|, IJJ, 1311) and close to the range of values {Cn/Ct « 9 - 14) 



found for nematodes crawling on substrates of varying stiffnesses (|2d ). 

The individual components of the drag force can be calculated by eval- 
uating the expressions enclosed by the curly brackets in Eqs. [11] and [T2j 
These expressions are a function of the normalized span length X = x/ Ryj 
and the normalized gap width e = ho/C, where C = {1/2){R~^ — R~^)R'^. 
Here, R^ = 40//m and Ra = 100/^m are the nematode's body radius and the 
groove's effective radius, respectively. The maximum value of span length 
Xm is 14 /im and ho = 0.85 fim. Hence, the values of the normal and tan- 
gential drag coefficients yield Cn = 222.0 and Ct = 22.1. To the best of our 
knowledge, the only other estimates available for Cn and Ct lie in the range 
ofCn^Ct^5- 40, with Cn/Ct < 2 (0). 



The effects of anisotropic drag forces on C. elegans' kinematic 
efficiency 

Traditionally, two main motility gaits (swimming and crawling) are reported 
depending on whether C. elegans is fully immersed in a liquid or moving on 
top of a surface. Parameters used to distinguish between these two motility 
forms have included metrics such as bending frequency / = 1 /T, amplitude 
A, and wavelength A, and values of suc h q uantities are found to be much 
lower for crawling than for swimming (jfil.bd). Our kinematic data also show 
similar trends (Table 1), and we find that the kinematic efficiency defined 
as U /c is eight times larger for crawling than for swimming. 

There are, however, some interesting similarities between crawling and 
swimming nematodes. In particular, the shape parameter value defined as 
\/A is nearly identical for both crawling and swimming gaits (Table 1). 
This constant value of \/A is likely to result from biological constraints on 
the muscle strain (jiol ). Since C//c is a function of Cn/Ct and \/A only (see 
Eq. [T3|) , we infer that the kinematic differences observed between these two 
motility gaits may be attributed to the different values of the ratio of drag 
coefficients experienced by C. elegans since \/A remains nearly constant. 

Overall, we find that anisotropic drag forces infiuence the undulatory 
locomotion in two aspects, i.e. the energy cost (Eqs. [TT]and[T2|) and the un- 
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Table 1: Kinematics for swimming and crawling nematode. 





Swimming 


Crawling 


Speed U (mm/s) 


0.36 ± 0.01 


0.26 ± 0.01 


Amplitude A (mm) 


0.26 ± 0.01 


0.10 ± 0.01 


Wavelength A (mm) 


2.12 ± 0.08 


0.79 ± 0.02 


wavespeed c (mm/s) 


4.20 ± 0.17 


0.38 ± 0.02 


Shape parameter (X/A) 


8.14± 0.29 


7.84 ± 0.33 


Efficiency (U/c) 


O.OSib 0.01 


0.67 ± 0.01 



dulation efficiency (Eqs. [T3]). The tangential drag coefficient Ct determines 
the nematode's required propulsion (mechanical) force and the energy cost 
because it is responsible for most of the translational resistance to forward 
motion, while the drag coefficient ratio C„/Ct determines the ratio of the 
propulsion energy dissipated along both the transversal and longitudinal 
directions (jlSl ). 



Calculating the nematode's internal bending force 

Bending is required for C. elegans to generate forward thrust. In this section, 
we calculate the total force necessary to bend the nematode's body using 
experimental data and knowledge of the drag forces. We compare the cal- 
culated nematode bending force patterns between crawling and swimming 
nematodes. A force balance on the nematode body shows that the bending 
force needed to overcome resistance arises from two components: one is the 
elastic resistance due to body bending, while the other is the viscous force 
distributed along the nematode body. In the limit of small undulating am- 
plitudes, the total muscle bending force F at position s = Sp can be obtained 
by summing up the elastic and distributed viscous drag forces acting along 
the nematode body, 

F{sp,t) = Kb^^^^^^n{sp,t) + fi J^' CnMs,t) + Ctut(.s,t) ds, (14) 

where s = and s = L represent the head and tail of the nematode body, 
respectively. The unit normal vector of the local body segment is given by 
n{sp,t). The first term on the right-hand side represents the elastic force 
required to bend the nematode body, while the second term represents the 
viscous force due to the liquid film. Altogether, internal muscle actuation 
must overcome both the elastic and viscous drag forces to sustain undulatory 
locomotion (see (js^ ) for details). 
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The total muscle force F is decomposed into normal and tangential di- 
rections. The tangential component of the force describes the resistance to 
changes in body length. The normal component of F(sp,t) describes the 
bending force, and it is given by 

Fb{sp,t) = F{sp,t)-n{sp,t). (15) 

The magnitude of the internal bending force Fh{s,t) can be determined 
at any time t and location s along the nematode body by combining exper- 
imentally measured nematode speed and bending curvature with individual 
values of the normal and tangential drag coefficients. To compute the in- 
ternal elastic contribution, we also need values of the nematode's bending 
modulus Ki), a material property that is typically difficult to measure. Here, 
we will use two previously estimated values of Kh, namely Kh = 4.2 x 10~^^ 
N m2 (30) and Kb = 9.5 x 10^^^ N m^ (@). The differences in the magnitude 
of Kb reflect the different assumptions in determining the spatial and tem- 
poral form of the nematode's (internal) active bending moment. Figure IHa, 
b) shows the magnitude of Fb using Kb = 4.2 X 10-^^ N m^ (30) for both 
crawling and swimming nematodes for three bending cycles. Figure HJc, d) 
shows the corresponding magnitude of Fb using Kb = 9.5 X 10"^^ N m^ (0). 
The solid black lines in Fig. U] correspond to regions where the nematode's 
curvature k is zero, and arrows show the direction of the bending wave. 

The bending force patterns presented in Fig. [J] are time-periodic and in 
general, spatially complex. For swimming nematodes with Kb = 4.2 x 10^^^ 
N m^ (l30l ). the bending force is larger near the head, while most of the 
bending force is accumulated around the mid-section for crawling C. elegans. 
For both swimming and crawling C. elegans, spatial patterns for the bending 
force are very similar to previously described gait-speciflc features of calcium 
signals in muscle (0). In particular, we note that the larger value for Kb 
(Fig. Iltjd) yields larger resistance to nematode bending. As a result, the 
spatial distribution of Fb is not identical for the two cases. 

Since the bending force patterns are time-periodic, the data shown Fig.[31^a,b) 
can be phase- averaged. Here, we focus our analysis on the Kb = 4.2 x 10~^^ 
N m^ ()30l ) case only. Figure [5]^a) and (c) shows the phase- averaged spatial 
distribution of Fb for crawling and swimming nematodes, respectively. Note 
that the shaded areas in Fig. [5] correspond to standard error of the mean 
(SEM). In general, crawling nematodes generate bending forces Fb that are 
fourfold larger than nematodes swimming in an aqueous solution. The peak 
values of internal bending force for crawling and swimming nematodes are 
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approximately 8.5 nN and 2.0 nN, respectively. As noted before, the bend- 
ing force during crawling appears concentrated in the middle section of the 
body. 

The phase- averaged bending force signal is shown as a function of one 
bending period in Fig. El^b) and (d) for crawling and swimming nematodes. 
The bending force signals are compared to the curvature signals, and mea- 
surements are made at s/L = 0.55 for both swimming and crawling. Results 
are shown for Ki, = 4.2 x 10^^^ N m^ only. We find that, for crawling, the 
bending force and bending curvature signals are not in phase with each other 
(Fig. [5)3). On the other hand, the force and curvature signals are nearly in 
phase for swimming nematodes. One possible explanation for the observed 
phase lag in crawling nematodes is the (much) larger viscous drag forces 
experienced by nematodes moving on wet surfaces. In the case where little 
viscous drag is present, the internal bending forces remain nearly in phase 
with the bending curvatures, as indicated by Eq. 14. In other words, a phase 
differences between the bending force and curvature is a consequence of en- 
hanced viscous drag forces; the larger the drag forces, the larger the phase 
difference between the internal bending forces and the bending curvatures. 
Finally, we note that for both gaits, the bending curvature signals are rather 
periodic, while the bending force signals exhibit unequal strength along ven- 
tral and dorsal sides. This unequal ventral/dorsal bending is most likely the 
result of an asymmetry of the ventral/dorsal neurological structures (jssi ) 
(see Fig. 4). 

We can further study the biomechanics of undulatory locomotion by com- 
puting the power dissipated by the nematodes during crawling and swim- 
ming. As shown above, the internal bending force is composed by the elastic 
and viscous forces, and the mechanical power dissipated by C. elegans during 
locomotion can be computed as 

P{t) = fi [ Cnulis, t) + Ctu^is, t) ds. (16) 
JO 

Since the elastic energy is conserved, energy dissipation is then due to the 
viscous forces. In Fig. [H we show the dissipated power for both crawling 
and swimming nematodes over one beating cycle, respectively. We find that 
the crawling nematodes dissipate more energy than swimming nematodes. 
The time-averaged power is 6.4 pW and 4.5 pW for crawling and swimming 
nematodes, respectively. Although the magnitudes in power are different, 
both the swimming and crawling gaits display similar temporal patterns. 
The phase differences between the power and the bending curvature tempo- 
ral signals are also quite similar for both motility gaits. 



11 



Discussion 



Uncoordination (unc) is a very common phenotype among genetic mutations 
in C. elegans. The affected genes may be expressed either in neurons, mus- 
cles or even in other tissues affecting the structures required for coordinated 
movement (jlTI ). For more than 30 years, different unc phenotypes of C. 
elegans were mostly observed on NGM (agar) plates. Here, we have quanti- 
fied the biomechanical properties of crawling C. elegans crawling using both 
experiment and modeling in an effort to provide quantitative tools for motil- 
ity phenotyping of mutant strains on wet agar gels. The analysis presented 
here focuses on wild-type C. elegans, although the methods can be easily 
extended to different nematode's strains and other organisms. We coupled 
kinematic data, such as nematode's speed and body curvature, with a hy- 
drodynamic model in order to obtain estimates of the (surface) drag forces 
acting on the C. elegans. Values of the surface drag forces were then used 
to calculate the nematode's internal (muscle) bending force, which are com- 
pared to recent calcium imaging data (26). Results on crawling nematodes 
are also compared to C. elegans swimming in water-like, buffer solutions 
in order to gain insight into the nematode's adaptive motility kinematics 
across different environments. 

Results show that the expressions for the surface drag forces along the 
normal and tangential directions of the nematode's body (Eqs. Illll2p share 
a similar form with drag forces obtained for swimming using resistive force 
theory (RFT) . Equations [11] and [12] indicate that resistance to locomotion 
can be quantified by the tangential drag coefficient Cj while propulsion effi- 
ciency is a strong function of the drag coefficient ratio Cn/Ct- We find that 
the drag coefficient ratio Cn/Ct is approximately 10 for nematodes moving 
on wet agar plates, in agreement with a recent study (0). Other investi- 
gations estimate values of Cn/Ct that range from less than 2 (27) to larger 
than 20 (0) for crawling nematodes. For comparison, the values of Cn/Ct 
range from 1.4 to 2 for swimming nematodes. By considering a thin liquid 
film between the nematode's body and the agar surface, we are also able to 
estimate the individual values of the normal and tangential components of 
the drag forces for crawling nematodes. Using Eqs. [TT] and [12] along with 
kinematic data, we find that C„ = 222.0 and Ct = 22.1 (Fig. [3]). For com- 
parison, we note that a recent theoretical analysis estimated C„ and Ct to 
be of 0(10), with a drag coefficient ratio less than 2 (j27l ). which is closer to 
swimming. 

In addition, our analysis of the drag forces suggests that gait modulation 
(crawling versus swimming) is probably a response to the external environ- 



12 



ment, or more specifically, to the drag coefficient ratio (Cn/Ct). That is, C. 
elegans seems to adjust not only the magnitude of the internal bending force, 
but more importantly its distribution along the body, to accommodate for 
changes in Cn/Ct and maintain its propulsive thrust, as shown in Fig. HI A 
recent investigation (jol) has shown that the power required for a nematode 
to swim in highly viscous fluids is orders of magnitude larger (5 nW) than 
found here for crawling, yet the ratio Cn/Ct for swimming remains constant 
at approximately 2 and is independent of fluid viscosity. Hence, nematodes 
maintain their distinct swimming gait even if the forces required to locomote 
are much larger than during crawling. Only when the drag coefficient ratio 
Cn/Ct is modified is that nematodes exhibit crawling-like gaits or behavior. 
Our present results suggest that C. elegans^ motility gait can be tuned and 
modulated by modifying the anisotropy of the surface drag coefficients. 

Knowledge of the individual components of the drag forces allows for 
calculating the internal bending force of both crawling and swimming ne- 
matodes (Fig. HI). We find that the bending force patterns are time-periodic 
and spatially complex. These patterns provide additional insight into the 
different motility gaits observed in C. elegans, namely swimming and crawl- 
ing. The temporal patterns of the bending force, for example, reveal a phase 
lag with respect to the bending curvature. This phase lag is a measure of 
the viscous drag force contribution, which is larger for crawling than for 
nematodes swimming in water- like buffer solutions. The spatial patterns 
show that large values of the bending force are concentrated near the head 
for swimming nematodes, while, for crawling, much of the bending force is 
concentrated in the nematode's mid-section. These calculated spatial pat- 
terns correlate rather well with regions of high calcium activity (j2^). The 
analysis and data presented here strongly suggest that while the motility 
patterns associated with crawling and swimming can be distinguished by 
distinct kinematics (Fig. [3] and Table 1) and bending force patterns (Fig. S]), 
changes in the environment's (external) drag coefficient ratio can also con- 
tribute to modulations in motility gaits (Eqs. 11 and 12). Further, the 
analysis is able to capture main biomechanical features of undulatory crawl- 
ing on lubricated surfaces (agar plates), and it can provide direct estimates 
of force and power expenditures that are useful for motility phenotyping 
activities such as in genetic and drug screening applications involving C. 
elegans. 
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Figure 1: (Color online) (a) The nematode C. elegans crawls on an agar 
plate. The red line shows the body centerline ('skeleton'). The black arrows 
show the unit normal and tangential vectors n and t along the nematode 
body. As the nematode proceed, it left over a groove on the surface, (b) 
Schematic of nematode cross-section during crawling. The fluid pressure on 
each side is p = pi and p = po, where po is the ambient atmosphere pressure. 
Other parameters are defined in the text, (c) The groove shape is assumed 
as a circular arc with radius Ra- The radius of the nematode is Rw Dashed 
lines are the parabolic curves and ha approximating the nematode and 
groove surface shape, respectively. 
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Figure 2: (Color online) The kinematics of C. elegans crawling on an agar 
plate. Body centerlines (or 'skeletons') of crawling nematodes in (a) labo- 
ratory frame and (b) nematode's frame of reference. The centroids of the 
skeleton in (a) are color-coded, and tracked to obtain the crawling speed U 
(see SV2). (c) Spatio-temporal contour plot of the nematode's bending cur- 
vature K, which display a traveling wave from head to tail, (d) Fast-Fourier 
transform of the curvature plot at varied locations of the nematode's body. 
The nematode's bending frequency is approximately 0.5 Hz and independent 
of position along the body. 
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Figure 3: (Color Online) The drag coefficients and the drag coefficient ratio 
(inset) as a function of the gap width ho. The blue (upper) and red (lower) 
lines are the normal and tangential drag coefficients, respectively. Inset: The 
normal and tangential drag coefficient ratio Cn/Ct as a function of Hq. The 
dashed lines show the estimated drag coefficient ratio Cn/Ct ~ 9.4 which 
corresponds to /iq ~ O.Sbfim. 
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Figure 4: (Color online) The normalized spatio-temporal bending force pat- 
terns for crawling (left-panel) and swimming (right-panel) C. elegans. The 
bending force calculated for (a,b) Kb = 4.19 x 10~^^ N m^ (|3ol ) and (c, d) 
Kb = 9.5 X 10"^'* N m^ (0). The overlaid black solid lines corresponds to 
regions where the curvature «; = 0. The dashed arrow shows the direction 
of the traveling wave. 
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Figure 5: Spatial (a,c) and temporal (b,d) patterns of the bending forces 
Ffy for crawling (top-panel) and swimming (bottom-panel) C. elegans for 
Kh = 4.19 X 10^^^ N m^. The shaded areas correspond to the standard error 
of the mean (SEM) (a,c) The spatial patterns are averaged over five periods. 
Most of the bending force is located near the head for swimming nematodes 
and around the body's mid-section for crawling. (b,d) Temporal patterns are 
obtained at position s/L = 0.55 for both crawling and swimming nematodes, 
and the dashed line corresponds to the bending curvature. 
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Figure 6: (Color online) Mechanical power dissipated during (a) crawling 
and (b) swimming. The solid red and dashed black curves are the dissipated 
mechanical power and bending curvatures. The shaded areas shown the 
standard error of the mean (SEM). Nematodes dissipate more energy during 
crawling compared to swimming. 
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